Condon Domain Phase Diagram and Hysteresis Size for Beryllium 
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The Condon domain phase diagram for berylhum is determined in magnetic fields up to 10 T 
and at temperatures down to 1.3 K using a standard ac pick-up coil method to measure the de 
Haas-van Alphen (dHvA) effect. The detection of the transition point from the homogeneous state 
to the Condon domain state (CDS) is based on the extremely non-linear response to the modulation 
field resulting from a small irreversibility in the dHvA magnetization. The experimental results are 
compared with theoretical predictions calculated from the Fermi surface (FS) of beryllium. The 
width km of the hysteresis loop in the CDS is measured in a wide temperature and field region. A 
model for the hysteresis size is proposed and numerically calculated for the whole phase diagram. 

PACS numbers: 75.45. -fj, 71.70.Di, 75.60.-d 
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I. INTRODUCTION 

CondonlJ predicted the formation of diamagnetic do- 
mains in non-magnetic metals, now known as Condon 
domains. A thermodynamic instability arises according 
to the Pippard-Shoenberg concept of magnetic interac- 
tiorPEl in the de Haas-van Alphen (dHvA) effect when 
the amplitude of the oscillatory magnetization signal be- 
comes large enough, i.e., the susceptibility 
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where M is the magnetization and B the induction. In 
this case the stability condition ^odH/dB = 1 — x > is 
not fulfilled for a certain interval of the applied magnetic 
field H in each dHvA cycle. For an infinitely long rod-like 
sample (demagnetizing factor n = 0) the system avoids 
the instability region by a discontinuous change of the 
induction B between the two stable states Bi and B2 
at a certain H = He- Both stable states have the same 
free energy (see also Fig. [9]) and the interval (_Bi,i?2), 
containing the instability, is forbidden. 

For a plate-like sample oriented normal to H {n — 1) 
the boundary condition B — (J-qH is required even within 
the interval Bi < fioH < B2, so that the induction B can 
not change discontinuously and the state with homoge- 
neous magnetization is impossible. The plate breaks up 
into regions of different magnetization with the induc- 
tions Bi and B2 ■ The volume fractions of the domains is 
adjusted in a way that the average induction B = ^qH 
is fulfilled for the whole sample^. For a sample with in- 
termediate demagnetizing factor < n < 1 the magnetic 
field interval Bi < fioH < B2 with domains decreases 
proportionally to n. Even for samples of arbitrary shape 
there is without doubt a non-uniform Condon domain 
state (CDS) with the same dia- and paramagnetic phases 



Bi and B2- However, the domain configuration is cer- 
tainly more complex. 

Up to now Condon domains have been observed by 
different experimental methods; by NMR^, /iSR spec- 
troscopj*^*^ and they were recently directly observed by 
Hall probes^. All these experiments have in common that 
two distinct inductions Bi and B2 or an induction split- 
ting AB = B2 — Bi are measured at a given applied field 
and temperature. 

Equation [T] defines the phase boundary between the 
uniform and the CDS which can be calculated using for 
example the Lifshitz-Kosevich (LK) formula for the os- 
cillatory dHvA magnetization resulting from the Lan- 
dau quantization of the conduction electrons in a metal. 
Theoretical calculations of this boundary exist for sev- 
eral metals^E]. However, the above cited measurements 
yielded only a few points in the (iJ, T) diagram where 
Condon domains were actually observed without a com- 
plete determination of the Condon domain phase dia- 
gram. An experimental determination of the CDS phase 
boundary, i.e., where AB approaches zero, is difficult and 
time-consuming^*^. Without doubt, another method for 
the experimental determination of the phase boundary 
is needed to obtain sufficient data for a comparison with 
the theoretical predictions. 

Recently, hysteresis was observed in the dHvA effect 
under the conditions of the CD9^. Due to the irre- 
versible magnetization, an extremely nonlinear response 
to a small modulation field arises in standard ac suscep- 
tibility measurements upon entering the Condon domain 
state. The out-of-phase part and the third harmonic of 
the pickup voltage rise steeply at the transition point to 
the CDS. Moreover, it was shown that the point {H,T) 
where the hysteresis arises is independent of the sample 
shape. The threshold character of these quantities allows 
to measure a Condon domain phase diagram with high 



precision and in a wide temperature and field range. This 
offers the possibihty for a more detailed comparison with 
the theoretical calculations. 

The CDS boundary of silver has been successfully de- 
termined with this method^^. The FS parameters of the 
nearly spherical FS of silver are well known. Therefore, 
for silver the CDS phase diagram in the {H, T) plane 
can be precisely predicted using the LK-formula with the 
Dingle temperature as a parameteil^. Good agreement 
was found with experimental data of the dHvA oscilla- 
tion amplitude in the homogeneous state^SJ and for the 
resulting CDS phase diagrarrP^. This demonstrated that 
the method using the nonlinear response for the deter- 
mination of the CDS phase boundary is correct. 

For beryllium the FS under c onsid eration consists of 
the well known electron " cigars""^^^. The curvature of 
the FS at the extremal cross sections is very small giving 
rise to a high dHvA amplitude. In addition there are two 
close dHvA frequencies of 940 T and 972 T which lead 
to a beat in the dHvA amplitude. Due to this frequency 
beat the CDS phase diagram is more complex compared 
to silver. Several models have been proposed to calcu- 
late a CDS phase diagram for beryllium; a 3-dimensional 
electron gas model using the LK-formula and a purely 
2-dimensional electron gas modeP. However, the calcu- 
lations were in contradiction to experimental data ob- 
tained by /xSR. This disagreement required new phase 
diagram calculations with a modified LK-formula for the 
intermediate, between 2- and 3-dimensional, FS of beryl- 
lium taking into account the real shape of the electron 
" cigars'^ii^. This model is in good agreement with at that 
time available /xSR data. Most recently another theo- 
retical calculation for the phase diagram was proposed 
using a different model representation of the quasi 2- 
dimensional Fermi surface of berylliurrP^. However, the 
very few experimental data available do not allow for a 
complete test of the recent calculations over the whole 
phase diagrarrPJSJ. 

In this work we determine the experimental Condon 
domain phase diagram for beryllium in the whole (-ff , T) 
plane for T > 1.3 K using the appearance of nonlinear re- 
sponse to an ac modu lation field for the detection of the 
phase boundarjff^H]^ Moreover, the width of the hys- 
teresis loop in the dHvA effect is measured as function 
of temperature and magnetic field in the CDS. Finally, 
a model for the origin of the hysteresis is proposed and 
numerically derived. 
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FIG. 1: (a) Pickup voltage normalized on the modulation 
level for high (400 /iT) and low (40 ^T) modulation level 
at 1.3 K. Due to hysteresis in the CDS the pickup voltage 
decreases if the modulation level is of the order or smaller 
than the hysteresis loop width, (b) Envelope of the difference 
between both curves of Fig. flta) showing the field regions 
where Condon domains exist for non-zero difference. The 
inset shows an expanded view of the difference signal with 
the appearance of hysteresis in a part of each dHvA period. 



The experiments were carried out at temperatures down 
to T = 1.3 K in a 10 T superconducting coil with a homo- 
geneity of better than 10~^ in a sphere with 1 cm diam- 
eter. Some experiments were made in a 16 T coil with a 
variable temperature insert to measure temperature de- 
pendencies at constant magnetic field. The modulation 
frequency was about 160 Hz. 



II. EXPERIMENT 

A standard pickup coil system was used for the ac 
measurements of the magnetic susceptibility. The results 
shown here were measured on the same rod-like sample 
as in Ref. [H] of sizes 8x2x1 mm^ with the long side 
being parallel to [0001]. The magnetic field is applied 
parallel to the long side of the sample. We found from 
our measurements a Dingle temperature of Tq = 2.0 K. 



III. RESULTS 

The phase transition point to the CDS can be deter- 
mined by several methods^ which are all based on the 
appearance of hysteresis in the dHvA effeclj^. Figure [l] 
shows the pickup voltage normalized on the modulation 
level for low and high modulation amplitude in a large 
magnetic field range at 1.3 K. Due to the hysteresis in 
the CDS the response to an ac modulation becomes ex- 
tremely nonlinear and the first harmonic amplitude of the 
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FIG. 2: Temperature dependence of the imaginary part of 
the pickup voltage for several modulation amplitudes at the 
paramagnetic part of a dHvA oscillation at the beat antinode 
at 3.6 T. The dashed line indicates the critical temperature 
Tc = 3.0 K where the Condon domain phase is entered when 
lowering the temperature. 




FIG. 3: Field dependence of the imaginary part of the pickup 
voltage measured at 2.0 K and 2.5 K with 40 /iT modula- 
tion amplitude. The steeply increasing negative amplitude is 
caused by the hysteresis in the CDS. The vertical dotted lines 
indicate the transition field to the CDS for 2.0 and 2.5 K. 



pickup voltage normalized on the modulation amplitude, 
usually corresponding to the susceptibility Xj decreases 
strongly at the paramagnetic part (x > 0) of every dHvA 
period. The amplitude damping is observed if the mod- 
ulation level is of the order or smaller than the width of 
the hysteresis loop. In absence of Condon domains the 
normalized pickup voltage is independent of modulation 
level. Therefore, the substraction of two curves, one mea- 
sured with high and the other with low modulation level, 
reveals the magnetic field ranges where Condon domains 
exist. In other words domains exist if the difference is 
greater than zero. Figure llTb) shows the envelope of this 
function at T = 1.3 K. The inset presents the detailed 
difference of both curves in a small region, which corre- 
sponds to a cut of the CDS phase diagram at T = 1.3 K. 
We see in Fig. W[h) that there is no difference between the 
normalized pickup voltages for magnetic fields exceeding 
6 T which implies that Condon domains disappear for 
fields higher than 6 T at 1.3 K. 

For the above described method to determine the CDS 
phase diagram, two field sweeps must be measured for 
each temperature. In order to detect even very small 
hysteresis the low modulation level must be as small as 
possible. Therefore, it is difficult to detect by this method 
the existence of Condon-domains at field regions where 
the hysteresis loop width is small. 

It was shown that the out-of-phase part and the third 
harmonic of the pickup voltage appear with threshold 
character whenever there is small hysteresis in the dHvA 
magnetizatiorpJJ. A measurement of one of these quan- 
tities offers therefore a simple alternative way to deter- 
mine with high precision the phase boundary of the CDS. 
Figure [2] shows temperature dependencies of the out-of- 



phase part of the pickup voltage for a large modulation 
level range at the beat antinode of the dHvA oscillations 
at a maximum of x near 3.6 T. At the critical temper- 
ature Tc = 3.0 K the out-of-phase signal drops down 
rapidly upon lowering the temperature. This indicates 
a sudden phase shift of the ac response with respect to 
the modulation signal. The phase shift is caused by the 
emerging hysteresis in the dHvA magnetizatiorPH_ y^Q 
find the same Tc, indicated by the dashed line in Fig. [2] 
for all modulation levels showing that the determination 
of the CDS phase boundary is independent of the used 
modulation level. 

In the following we will use a modulation amplitude of 
40 ^T. This value is sufficiently small to insure that the 
dHvA period is always much bigger than the modulation 
amplitude h even at low magnetic fields {&t jIqH — 1 T 
the dHvA period of beryllium is about 1 mT). On the 
other hand the ac response is still easily detectable. We 
note that if the modulation amplitude is of the order of 
the dHvA period the imaginary part and the third har- 
monic of the pick-up signal show up even in the absence 
of hysteresis-. 

Figure |3] shows magnetic field dependencies of the 
imaginary part of the pickup voltage measured at 2.0 K 
and 2.5 K. The phase of the lock-in amplifier is adjusted 
such that the signal due to the sample susceptibility is 
mainly in-phase. dHvA oscillations of small amplitude 
similar to the waveform in the inset of Fig. [T] arc visible 
in the out-of-phase signal with an amplitude bigger than 
previously observed in Ref . [HI due to the increased mod- 
ulation frequency of 160 Hz compared to 21 Hz in Ref.fTTl 
The higher eddy currents explain the appearance of an 
out-of-phase signal for an homogeneous magnetization. 
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FIG. 4: Temperature dependence of the imaginary part (a) 
and the third harmonic (b) of the pickup voltage measured in 
the paramagnetic part of a dHvA oscillation at the beat node 
for 2.5 T with 40 jj,T modulation amplitude. 



However, Fig. [3] shows threshold character in the aris- 
ing of negative peaks at magnetic fields where hysteresis 
occurs indicating the transition to the CDS around the 
negative peaks in each dHvA period. We see in Fig. [3] 
that the negative peaks appear at 2.0 K at lower fields 
than at 2.5 K. At temperatures above 3.0 K all negative 
peaks disappear and only the small dHvA oscillations due 
to eddy-current effects remain. 

The amplitude of the negative peaks depends on the 
modulation level, the hysteresis loop width at the partic- 
ular magnetic field, and on the amplitude of the in-phase 
part of the pickup voltage, i.e. the susceptibility. Even 
though the peak amplitude seems to be correlated with 
the Condon domain phase diagram being stronger fur- 
ther away from the phase-diagram boundary, we extract 
from these data only the magnetic field values for which 
the negative peaks appear for each temperature in order 
to construct the phase diagram in the next section. The 
negative peaks arise with threshold character and Fig. [3] 
shows that the CDS phase boundary can be determined 
with a precision of about one dHvA period. 

It was reported in fiSR studiea^E^El that Condon do- 
mains occur also at the beat nodes of the dHvA oscil- 
lations around 2.0 T and 2.7 T for 0.5 K and 0.8 K, re- 
spectively. However, there are only a few temperature de- 
pendencies of the induction splitting available from /iSR 
measurements. In other words, the reported tempera- 
tures do not represent necessarily the CDS phase bound- 
ary for these fields. Figure |4] shows the temperature de- 
pendence of the out-of-phase part (a) and the third har- 
monic (b) of the pickup voltage at the beat node at 2.5 T. 
A sharp transition at 1.5 K is visible in both traces which 
indicates that hysteresis arises at this temperature. This 
means that Condon domains appear indeed at this beat 
node and the CDS phase boundary is at 1.5 K for 2.5 T. 
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FIG. 5: Phase diagram in the {H, T) plane for beryllium. 
Scatter points indicate the position of the phase boundary at 
beat nodes and antinodes determined by temperature sweeps 
(circles and triangles) like in Figs. p] and |4] and by field sweeps 
(squares) like in Fig. ^ The solid lines are guiding lines of the 
phase boundary for beat nodes and antinodes, respectively. 
The dashed lines show for comparison the recent theoretical 
calculation of Ref. I6I15I The dotted lines indicate the enve- 
lope of the beating substructure of the phase diagram shown 
in detail in Fig |6] 



IV. PHASE DIAGRAM 

We have seen that the CDS phase boundary can be 
determined with high precision using nonlinear response 
measurements. Due to the hysteresis the out-of-phase 
signal of the pickup voltage drops sharply. This was 
measured either at a fixed magnetic field like in Fig. [2] 
and |4] or at fixed temperature as function of magnetic 
field like in Fig. [3] All data is compiled to obtain a 
complete Condon domain phase diagram in Fig. [51 The 
solid lines in Fig. [s] are extrapolated guiding lines to the 
{H, T)-values obtained for the beat antinodes and nodes, 
respectively. These lines are the envelopes of a substruc- 
ture consisting of a beating pattern of sharp needle-like 
domain regions as shown in Fig. [6] where the inset reveals 
the Condon-domain regions in two successive dHvA pe- 
riods. For magnetic fields in between these needle-like 
regions the sample is in the homogeneous state. We see 
in Fig. [6] that Condon domains appear first for magnetic 
fields around a beat antinode where the dHvA ampli- 
tude is higher. When cooling down the CDS field range 
extends gradually around the antinodes. 

The condition (see Eq. 1) that a Condon domain state 
occurs in a dHvA period is independent of the demagne- 
tization factor. We have found in a test on a plate-like 
sample with the same Dingle temperature that the ob- 
tained phase diagram is indeed independent of the sample 
shape. However, this only holds for the envelope of the 
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FIG. 6: Expanded view of Fig. |5] around the beat maximum 
at 4.8 T showing the beating substructure of the phase dia- 
gram. The shght steps in the envelope of the height of the 
needle-hke stripes results from the limited number of mea- 
sured temperatures. The expanded view in the inset shows 
the detailed phase diagram in two successive dHvA periods 
obtained by field sweeps at constant T. The larger round 
circles are obtained by T-sweeps at constant H. 



phase diagram (solid lines in Fig. I5| . The substructure 
depends on the sample shape as the needle-like regions 
(inset of Fig. [6]) are much broader for a plate-like sample. 
The reason for this is that the field range within a dHvA 
period where domains arise scales with n and is therefore 
more extended in a plate-like sample^. In other words, 
the envelope of the phase diagram in the {H, T)-plane is 
independent of the demagnetization factor, but not the 
detailed field region within a single dHvA period. 

The p hase di agram in Fig. [5] agrees with all reported 
/iSR dateP^SElliZI^ j^ particular, the observed induction 
splitting disappeared at the beat antinode near 2.6 T for 
temperatures higher than 3.0 K. We examined the same 
beat maximum and found a critical temperature of 2.9 K 
for our sample. 

In Fig. [5] the experimental phase diagram is compared 
with the calculations made with the modified LK-formula 
in Ref. I6I15I The overall shape of the calculated antinode 
and node envelope curves is similar to the experimental 
result. However, there is a clear discrepancy between the 
predicted temperature and magnetic field ranges of the 
CDS and the ones we find experimentally. We observe 
an upper critical field of about 8 T (extrapolation of the 
guide to the eye for the antinodes in Fig. [s] to T = 0) 
above which domains disappear for all temperatures for 
our sample with To — 2.0 K. We see in addition that 
Condon domains continue to exist at higher temperatures 
down to lower fields compared to the phase-diagram cal- 
culation and that domains exist at the beat nodes up to 



higher temperatures than predicted. A reason for the 
discrepancy in the temperature-field values of the phase- 
diagram boundary can be possibly related to the strong 
magnetostriction effects in beryllium (see discussion be- 
low). 



V. HYSTERESIS LOOP SIZE 

Another interesting question is the dependence of the 
hysteresis loop size hm from temperature and magnetic 
field. At the CDS phase boundary hm vanishes. How- 
ever, the temperature and field dependence of hm in the 
CDS might give information about the nature of the ir- 
reversible magnetization. We suppose that the hysteresis 
itself and its size are mainly caused by irreversible domain 
wall motion or by nucleation of new domain fractions. In 
the following we will measure the temperature and field 
dependence within the phase diagram. 

As it was shown earliei'ii' hm can hardly be measured 
directly by Hall probes because of its small magnitude. 
However, hm can be indirectly determined by analyz- 
ing the response characteristic to an ac modulation field. 
As shown above, the normalized pickup voltage decreases 
strongly if the modulation amplitude decreases below the 
hysteresis size (see Fig. I7|. All measurements were made 
in the center of a dHvA period, i.e., in the center of 
a needle-like stripe of the phase diagram (see Fig. |6|. 
From Fig. [7] the order of magnitude of hm can be es- 
timated. However, it is not obvious which modulation 
level corresponds actually to the real hm which would 
be observed with Hall probes. A comparison with direct 
measurements by Hall probes under the same conditions 
yields good agreement if we chose for hm the onset of the 
decrease in the normalized pick-up voltage indicated by 
arrows in Fig. [7] All data points obtained in this way 
are presented in Fig. [8] We see a more or less linear tem- 
perature dependence of hm far enough from the phase 
boundary. Moreover, hm is practically independent of H 
at the lowest measuring temperature far from the phase 
boundary. 



VI. MODEL FOR HYSTERESIS LOOP SIZE 

Hysteresis in the CDS is certainly due to irreversible 
domain wall motion or rearrangement processes of the 
respective domain volume fractions upon field variation. 
In the following we analyze the shape and amplitude of 
the energy barrier in the domain wall between the phases 
with the inductions Bi and i?2. 

We can write the potential il as the sum of the contri- 
bution of the magneto-quantum oscillations (taking only 
the first harmonic of the LK-formula) and the magneto- 
static energy 
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FIG. 7: Normalized pickup voltage as function of modulation 
level for several temperatures in the center of the dHvA os- 
cillation at 2.6 T corresponding to the center of a needle-like 
stripe like shown in Fig. |6] The arrows indicate the chosen 
hysteresis sizes hm- 
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FIG. 8: Hysteresis amplitude hm as function of magnetic field 
and temperature (data points), determined like in Fig. IT] as a 
measure for the hysteresis width. The parabola indicates the 
region of the CDS phase diagram in the [H, T) plane (solid 
line of the antinodes in Fig. Isl. 



with the dHvA frequency F and the oscillation ampli- 
tude a{H, T) given by the LK-formulsP. For inductions 
B close to the magnetic field H we can develop Eq. [2] 
setting B — ^qH + b 
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FIG. 9; Schematic representation of the Lifshitz-Kosevich 
potential in the CDS showing two minima corresponding to 
the domain states with the inductions Bi = fioH + bi and 
B2 — noH + b2 and a potential barrier AE in the domain 
wall between the domains. 



a{H, T) is big enough, which corresponds to the condition 
in Eq. [l] then 51 has two minima at &i and 62. Figure [9] 
shows schematically the potential under this condition. 

The states between Bi and B2 have extra energy and 
are not stable. However, in the domain wall the induction 
B has to cross all values between Bi and i?2. Therefore, 
there is an energy barrier (see Figure l9| separating the 
states with these inductions whose amplitude /S.E can 
be calculated. We would expect that /S.E scales with 
the hysteresis width because this energy barrier must be 
overcome when the domain distribution changes under 
variation of the applied field. 

First, the inductions values Bi and B2 of the domains 
are found by minimization of the free energy 
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The induction difference A_B between the domains is cal- 
culated as a function of H and T in Fig.[TO] Here, and in 
the following calculations we use the 'cylinder' model of 
Ref. 15 which approximates the cigar-like Fermi surface 
shape of beryllium with a cylinder to determine x(ff, T) 
and a{H,T). This idealized model gives a reasonable 
upper limit for the phase boundary of the domains and 
should be sufficient to get an idea of the overall behavior 
of AE as a function of H and T. 

Once &i and &2 are known we can calculate the ampli- 
tude of the energy barrier AE using the above formulas 



with the dHvA period p = {^qH)^ /F. If the amplitude 
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FIG. 10: Calculation of the induction difference AS between 
the domains as a function of H and T for a Dingle temperature 
Td = 2.0 K. 



A simpler expression for AE can be given in good ap- 
proximation taking into account that the shape of ri(6) 
between bi and 62 is very similar to the cosine-function 
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Taking the second derivative of this function with respect 
to b and taking into account that the curvature of fl is 
1/fiQ — X at 6 = 0, (here x is positive and ^oX > Ij '^6 
find the following expression 
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This expression for AE is simi lar to the do main wall sur- 
face energies calculated earliepl^l^. Figure 
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shows the 

numerical calculation of AE as a function of temperature 
and magnetic field. 

We see a qualitative agreement between the calculated 
behavior of AE with the observed width of the hysteresis 
loop in Fig. |8] Both show a maximum at approximately 
the middle of the magnetic field range of the CDS phase 
boundary at lowest temperature. There is also good 
agreement with a recent theoretical calculation where the 
Rayleigh model was used for the hysteresia^SI. However, 
we notice that the energy barrier model shows a clear 
difference from our data near the phase boundary, h^ 
varies rather steeply when approaching the phase bound- 
ary whereas AE increases only gradually at the phase 
boundary. One can say that our model describes well 
the Condon domain hysteresis far from the phase bound- 
ary where the irreversibility results mainly from domain 
wall motion. We suggest therefore that close to the phase 
boundary a mayor contribution to the hysteresis is due to 
the nucleation of new phase fragments in tubular form. 
The region where the tubular structure undergoes a tran- 
sition to a laminar one with stripes could be very narrow 
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FIG. 11: Numerical calculation of the energy barrier AE be- 
tween the domains with the respective inductions B^ and B2 
as function of H and T for a Dingle temperature Td ~ 2.0 K. 
The solid line indicates the Condon domain phase boundary 
using the 'cylinder' model of Ref. 1151 



compared to the intermediate state of type-I supercon- 
ductorg^il. 



VII. DISCUSSION 

We can summarize two main differences in the behav- 
ior between beryllium and silver. First, as it was shown 
recently by using local Hall probes^, Condon domains do 
not emerge to the sample surface in beryllium, appear- 
ing only inside the bulk. For silver the measured induc- 
tions values of the respective domains are practically the 
same inside the bulk and on the surface indicating that 
domains emerge completely to the surface. Second, as 
it is shown here, the experimental CDS boundary lies 
inside the closest model and the discrepancy increases 
for higher magnetic fields unlike to silver where at least 
up to 30 T good agreement was observed with theoreti- 
cal calculationJi^l. We propose that the reason for both 
findings could lay in the fact that the dHvA effect is al- 
ways accompanied by magnetostriction oscillations. In 
the particular case of the CDS this means that domain 
formation gives rise to different, act ually opposite, defor- 
mations in the neighboring domainP^BU 

This deformation varies across the domain walls be- 
tween neighboring domains and requires extra elastic en- 
ergy in the domain walls and on the surface. Moreover, 
the magnetostriction amplitude increases with magnetic 
field and the amplitude is actually very big especially for 
beryllium. Beryllium has in comparison to silver a much 
higher Young modulus and the deformation under the 
dHvA effect is anisotropic. This idea can qualitatively 
explain the discrepancy in the phase diagram between 
theory and experiment at high magnetic fields. Neverthe- 



less, the discrepancy remains for low magnetic fields and 
for the nodes where the experimentally obtained bound- 
aries are above the calculated ones. It would be of in- 
terest to include the influence of magnetostriction on the 
calculated phase boundaries. 



VIII. CONCLUSION 

We have measured a complete Condon domain phase 
diagram for beryllium at temperatures down to 1.3 K 
and magnetic fields from 1 T up to 10 T. The method 
based on the detection of the nonlinear response to an 
ac modulation field provided also information about the 
substructure of the phase diagram due to the dHvA fre- 
quency beat in beryllium. The measurements agree with 
all data obtained by /iSR. Moreover, we have checked 
that the obtained phase diagram is independent on the 
sample shape. The method can be easily applied to sam- 
ples with other Dingle temperatures and other metals. 

The hysteresis loop size was measured in a wide region 



of the CDS phase diagram. In the middle of the phase 
diagram, far enough from the boundary, the hysteresis 
loop width increases linear with decreasing temperature 
and it is almost constant with magnetic field. Finally, a 
model for the origin of the hysteresis is proposed. The 
induction difference between the different domains is nu- 
merically derived and the height of the energy barrier 
separating these two states of induction is calculated. We 
found that the calculated energy barrier scales well with 
the observed hysteresis loop width besides in the region 
close to the phase boundary. We suggest that close to the 
phase boundary the domain wall motion is not the only 
reason for the observed hysteresis and that the process of 
filamentary nucleation of the newly created phase must 
be taken into account. 
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